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Abstract — Communication over a multiple access channel is 
considered. Each user modulates his signal as a superposition of 
redundant data streams where interconnection of data bits can be 
represented by means of a sparse graph. The receiver observes 
a signal resulting from the coupling of the sparse modulation 
graphs. Iterative interference cancellation decoding is analyzed. 
It is proved that spatial graph coupling allows to achieve the 
AWGN channel capacity with equal power transmissions. 

I. Introduction 

The effect of spatial graph coupling, first discovered for low- 
density error correction codes, attracts increasing interest from 
the communications research community. Iterative decoding 
of block low-density parity-check (LDPC) codes, despite its 
known strength, fails to achieve the same limits as the optimum 
maximum likelihood (ML) decoding. It has been shown, how- 
ever [1], and then rigorously proved [2], that iterative decoding 
thresholds of convolutional (or spatially coupled) LDPC codes 
coincide with the ML decoding thresholds. A coupled code 
is constructed by copying the protograph of an initial block 
code a number of times and then connecting the neighboring 
copies to form a chain. The graphs located at the ends of 
the chain contain variable nodes with fewer constraints, since 
these graphs are connected to their neighbors from only one 
side. As a result, iterative decoding progress initiates at the 
ends of the chain, due to the effect of the slight irregularity, 
and then propagates through the entire chain. The principle of 
spatial graph coupling has proven to be applicable to several 
other areas, including multiple user detection, compressive 
sensing, and quantum coding. 

In this paper, we focus our attention on spatial graph 
coupling in relation to multiple access communications. The 
capacity and the achievable rate region for the multiple access 
channel (MAC) are well known. The corner points of the rate 
region polytope can be achieved by successive (onion peeling) 
decoding [4] while the middle part is achieved by rate splitting 
and/or time sharing added to the onion peeling [5]. However, 
communication approaches that allow for more robust and 
less complex joint parallel detection/decoding, such as code- 
division multiple access (CDMA), have had limited success in 
achieving the inner points of the MAC rate region. The equal 
power user case is typically the most difficult, since all users 
happen to be operating under same conditions, and there is no 



structural irregularity to initiate the decoding convergence. It 
has been shown that regular random CDMA can only support 
a fixed system load 1 equal to a = 1.49 [8] with equal power 
users, while sparse synchronous CDMA can only support 
a = 2.07 [6]. 

In this work we consider a multiple access multiuser sit- 
uation where each user modulates his signal as a sum of 
redundant, independent, equal power, data streams, as de- 
scribed in [6] [7]. The modulated signal allows for sparse graph 
representation similar to the LDPC code's Tanner graph. The 
sparse graphs of the transmitted signals couple together "in the 
air," and the receiver observes a coupled sequence of graphs 
which is later processed by an iterative detector/decoder. We 
focus on a communication scenario in which all users transmit 
an equal number of equal power BPSK data streams that 
are encoded by error correction codes optimal for the binary 
input additive white Gaussian noise (AWGN) channel. We then 
prove that the (real-valued) AWGN channel capacity can be 
achieved exactly by such a system. 

II. System Model 

We consider the generalized modulation type format [7]. 
Each user modulates a number of independent data streams 
and transmits their sum. One data stream is modulated as 
follows. First the data is encoded by a binary error control 
encoder to produce a binary data stream {ui}f" =1 . The stream 
is first duplicated M times resulting in M identical data 
sub-streams, {u m j}^ =1 , m = 1,2, ... ,M, Each sub-stream 
{ u m,i}[Li i s permuted by an interleaver 7r m to produce 
{um,i}iLi> m G {I? 2, . . . , M }. At the next step each bit 
of {u m ,i}fLi is multiplied by an A*" -dimensional signature 
vector s m . Signature vector maps the data into the available 
resource space. This operation is similar to spreading in a 
CDMA system. Each of the M sub-streams is multiplied by 
an amplitude \J P/Al and they are all added up to produce 
one modulated signal stream 
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vi = ^ Y —u m ,is m 1 = 1,2, ,..,L. (1) 
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1 The system load a is defined as the ratio of the number of supported users 
K to the number of available signal dimensions N (or chips in CDMA), i.e, 
a = K/N. 
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Fig. 1. Transmitted packets coupling in the channel, W = 2. 



Each signature sequence s m — {s m ,i, s TO ,2) • • • i Sm,N} is 
real and energy normalized such that E|s m j\ 2 = ^, i = 
1,2, . . . , N. We also require 2 Es„ ll s'^ l — l/N, mi ^ m^. 

We consider the case in which 2W + 1 users transmit 
Kj (2W+1) data streams each, producing total system load of 
a = K/N. The power of each modulated stream is P = 1 and 
communication happens over the real-valued AWGN channel 
with noise power a 2 . We focus on a scenario in which the 
packets of the users are received asynchronously as shown 
in Fig. 1. Each packet is divided into 2W + 1 sections of 
length (2W + l)/L bits each. At time slot t = -W the first 
user starts to transmit his packet and transmits one section. 
At time t = —W + l the second user joins and transmits 
one section while the first user transmits his second section. 
At time t = —W + 2 the third user joins and so on. Once a 
user finishes transmitting his packet he immediately starts to 
transmit his next packet. 

We consider two types of receivers. The first receiver is a 
two-stage receiver where the first stage is an iterative parallel 
interference cancellation (PIC) which layers the received data 
streams and the second stage is an error control decoding per- 
formed for all layered streams independently and in parallel. 
The first, PIC stage, can be described as follows. At every 
iteration the received signals are filtered using the signature 
sequences Sk, m - Assuming a simple matched filter, the filtered 
signal of data bit ; from data stream k e {1, 2, . . . , K}, and 
sub-stream me {1,2,..., M} equals 

2,lSk,m ^ ^ ^fci,mi,^);i,mi) 
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where Ik, m .t is the interference, suppressed by the filtering. 
Each bit Uk,i has M replicas Uk. m ,i which are dispersed along 
the packet. The power x m of the interference term Ik.m.i de- 
pends on the section of the packet where replica Uk, m ,i belongs 
to. The signal power of each replica is 1 jM. We estimate the 

2 The following property can be easily guaranteed by choosing s m i £ 
< — ^7=, ^7^} Bernoulli with probability 1/2. However, other signature 
sequence choices are also possible. 



transmitted bits by summing up weighted received replicas 
1k,m,i and applying the conditional expectation estimate 
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The optimal weighting coefficients 



ensure that the signal-to-noise-and-interference ratio (SINR) of 
the weighted sum inside the parenthesis in (3) is maximized 
and equals 

The bit estimates (3) are used to reconstruct the entire trans- 
mitted signals (1) and perform interference cancellation by 
subtraction of the reconstructed signals from the received 
signal. The iterative process is repeated i times until the SINR 
Zi of the data bits at iteration i becomes greater than the 
decoding threshold 9 of the error control code used for their 
encoding. The error control decoding (second stage) can than 
be applied to all the packets (received at given time) for the 
final error correction. 

The second receiver, which we call a modified successive 
interference cancellation (SIC) receiver is working as follows. 
It starts with the first stage of the two-stage receiver described 
above. However, the error correction decoding of the data 
streams for which Zj > 6 is performed after each interfer- 
ence cancellation iteration. These decoded streams are then 
removed from the interference pool. The modified SIC receiver 
is slightly more complex since it requires a feedback loop 
between the decoder and the interference canceller. However, 
the modified SIC allows for simpler analysis. 

III. Analysis 

We start with an analysis of the iterative interference can- 
cellation, i.e., the first stage of the two-stage receiver. We track 
the evolution of the noise-and-interference variance which we 
denote by x\, and the SINR of the data bits which we denote 
by z\. The upper index t indicates the time slot of interest 
and i is the iteration number. A discrete system of equations 
describing the evolution of x\ and z\ have been derived in [3] 
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where the function <?(•) given by 
g(a)=E (l -tanh(a + ^ v / a))' 



i > 0,t > 1 (5) 



i > 0,t > W + 1 (6) 



£~JV(0,1) 



determines the mean squared error of the data bits with SINR 
a and Af(0, 1) denotes a standard normal random variable. 
Function g(-) is a continuous differentiable function that is 
strictly decreasing from g(0) = 1 to g(+oo) = 0. Variable 
z\ defines the SINR of the packet starting at time t — W 
and finishing at t + W, i.e., centered at time t. SINR z\ 
is computed from the SINRs of the M partitions uniformly 
distributed along the packet. The SINRs -4+r of the partitions 

x i-l 

are averaged in (6), just like in (4). Packets that start between 
time t — 2W and t contribute to the noise-and-interference 
variance at time t in (5). 

The iterative system is initialized by 



+oo, 
0, 



t < w 
t>w 
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(7) 



since the packets that start before time t = are completely 
known (or absent) and, therefore, their SINR is +oo. These 
packets are centered at time before W + 1. The packets 
centered at time t = W + 1 and after are completely unknown. 
Hence, their SINR is 0. 

Considering W — > oo and normalizing the packet length to 
1 the equations (5) and (6) can be transformed to 
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n / g(z t j +T )dT + a 2 
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dr 
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t>~,i>0 (8) 



t>0,i>0 (9) 



with a shift of the time origin from 1/2 to which is done 
for notation simplicity. The initialization is done using a step 
function 

+oo, t < 
0, t > . 



(10) 



The SINR of the data bits at iteration / of the interference 
cancellation stage equals z\. The second, error correction 
decoding stage, of the two-stage reception is completed suc- 
cessfully iff z\ > 9 where 9 is the error correction code 
threshold. 

Contrary to the two-stage receiver, the modified SIC receiver 
is performing decoding after every interference cancellation 
iteration. SIC can be described by (8) together with 
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replacing equation (9). We notice that the operator F(-) 
defining the evolution of the SINR throughout the iterations 
via (8) and (11), 

F{z\) , (12) 



is monotone, i.e., if Z\(t) < f° r t € (— oo, +oo), then 

F( z i{t)) < F(z 2 (t)) for t g (-oo,+oo). 



Theorem 1. Spectral efficiency of 
«Cbiawgn 
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is achievable by the modified SIC receiver for any a G [0, oo) 
and e > 0. 

Proof: We proceed with (10) and (8) calculating 
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x l - I aa (z T ) dr + a 2 — I at + ^ + <e hil] 
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that after the application of (9) implies 
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-S + o(S), t e [0,S] (14) 



a 2 (a + a 2 ) 

for small 6 > 0. Let us assume that the binary error correction 
codes for all data streams are of rate 

Cbiawgn f - In —^5 57 — ■ — ^r6 + o(S)) , (15) 

\a o cr^{a + a z ) ) 

i.e., these codes have threshold 
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In this case (11) implies 
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t < S 



6 + o(5) 



indicating that the packets centered at t G [0, S) are success- 
fully decoded and cancelled. The situation will repeat at the 
next iteration since Zg IC x > zff~ 5 , t G (— oo, +oo) implies 
z sic 2 = +°° f° r * e [^s 25) from the monotonicity of the 
SINR evolution operator (12). Thus, the convergence of the 
SINR to +oo will continue at the speed of at least 5 per 
iteration. The total achievable sum rate is (see (15)) 

I ' ( 1 . a + a 2 a /r , N 

a Cbiawgn - In = T — ■ — ^-o + o(d) 

By choosing appropriately small 5 we obtain the statement of 
the Theorem. ■ 
Let us define the total system SNR parameter s = The 
next theorem states that for any fixed s the AWGN channel 
capacity is achievable by the modified SIC decoder. Let us 
consider a 2 = a/s, i.e., we are keeping s fixed. We denote 
the limiting spectral efficiency of the system considered in 
Theorem 1 and the corresponding SNR per bit by 

Ces(a, s) = aC B iAWGN f — ln(l + s) ) (16) 
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We also define the real-valued AWGN channel capacity 
Cawgn (tTo) f° r gi ven jjg as the ro °t °f the equation 

C ™(t)^ l0g2 ( 1 + 2CAWGN (f)t) ' (18) 
Finally, we state the main result of the paper. 
Theorem 2. 



lim C e ff(a, s) = lim Cawgn ( —(a,s] 
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log 2 (l + s) 



Proof: The capacity of the binary input AWGN channel 
for SNRs 7 < 1 can be bounded as follows [9] 

^^ 7 + A l7 2 <C B iawgn(7) < ^^7 + ^27 2 (19) 

where A\ < A 2 are constants. Therefore, application of (19) 
to (16) implies 
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while application of (19) to (17) implies 

s Eb, 

lo g2 (l + S ) + 2A 2 ^-N ^ S) 

s 

< 



log 2 (l + s )+2A 1 ^sM 
Taking the limit in the above inequalities we obtain 

lim C e ff(a, s) = - log 2 (l + s) 

a—>oo Z 

,. Eb s 

lim — — (a, S) = : 7- r 

a^coiV ' log 2 (l + s) 

Equations 22 and 23 imply equality 
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which ensures (see (18)) 
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The Theorem is proved. ■ 
Spectral efficiency C e g(a, s) achieved by the modified SIC 
receiver (see Theorem 1) is plotted in Fig. 2. The three 
magenta curves correspond to a = 10, 100, and 500 (from 
bottom to top). For each curve a is kept constant while s 
is varying. The channel capacity Cawgn is given by the 
blue curve. Finally, the brown curve plots spectral efficiency 
achieved by the two-stage receiver. 




Fig. 2. Achievable spectral efficiency for the modified SIC receiver 
(magenta), two-stage receiver (brown) and the AWGN channel capacity (blue). 



IV. Conclusions 

In this paper, it has been proven that sparse graph modula- 
tion with spatial coupling can achieve AWGN channel capacity 
under modified SIC reception. The sparse graph modulation 
format is based on superposition of low-rate redundant data 
streams which can be easily designed to operate near binary 
input AWGN channel capacity. Numerical results also demon- 
strate that the two-stage reception where parallel interference 
cancellation is followed by parallel error control decoding is 
capable of operating within ldB of the channel capacity. 
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